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Abstract. Brown dwarfs bridge the gap between the stellar and planetary mass 
regimes. Evolving from conditions very similar to the lowest-mass stars, the atmo- 
spheres of older brown dwarfs closely resemble those expected in close-in extrasolar 
giant planets, and with cooler BDs still being discovered, more and more approach the 
properties of gas giants at wider separation. Interpreting the spectra of BDs is there- 
fore a crucial step towards understanding and predicting the spectral and thermal 
properties of EGPs. 

Essential properties of substellar atmospheres are massive molecular line-blanketing 
and the condensation of species with decreasing T e B , changing the chemical equilibrium 
composition and causing absorption from dust grains. More complex details involve 
the distribution of dust clouds over the surface giving rise to temporal variability, and 
possible deviations from chemical equilibrium conditions. In the case of close-in EGPs 
and some BDs in binary systems, the effect of irradiation from the primary significantly 
affects the spectral properties and thermal evolution. 

1 Introduction 

Observational efforts during the past decade have brought marked progress in 
characterising the lower end of the Main Sequence down to and beyond the 
hydrogen-burning limit, and identifying a class of unambiguously substellar ob- 
jects. Brown dwarfs (BDs) are commonly defined as compact objects with a mass 
below the minimum for sustaining equilibrium hydrogen fusion (~ 0.07Mq). 
Since the detection of the first bona-fide BDs |1I2I3| and the discovery of a 
Jupiter-mass companion to 51 Peg by Mayor & Queloz 0], the field of substellar 
astronomy has seen the direct detection of more than 300 ultracool dwarfs and 
the indirect detection of 110 planetary mass objects in orbit around other stars. 
While there is no agreement yet whether brown dwarfs should be distinguished 
from extrasolar giant planets (EGPs) based on their formation history, or by 
defining the minimum mass for deuterium burning (~ 13Mj) as the lower limit 
of the BD mass range, in terms of atmospheric properties and spectroscopic ap- 
pearance there is a smooth transition, and some overlap, between both classes 
of objects, regardless of definition. 
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The spectroscopic characterisation of these sources of 2200 K > T c g > 700 K 
has prompted the introduction of the two new spectral classes L and T. While the 
original description of the L sequence was based on optical spectral properties 
j5l(ij . the definition of a classification scheme for the latest L-type and the T 
dwarfs has required the use of near-IR spectra, owing to the fact that the lion's 
share of emitted flux and of spectroscopic characteristics in these objects is found 
at A > 0.95 nm [715] . 

While the youngest / most massive brown dwarfs jointly populate the early 
L (and late M) classes together with the least massive hydrogen-burning (very 
low mass-, or VLM-) stars, more evolved brown dwarfs of late-L type have effec- 
tive temperatures similar to those of the EGPs closest-in to their primary star 
("Ho t Jupit ers", e.g. 51 Peg, HD 209458 0). The coolest T dwarfs detected to 
date |10lll| exhibit atmospheric properties close to those anticipated in EGPs 
at larger separation, and thus are crucial test cases for our understanding of 
planetary atmospheres. 

The spectral appearance of substellar objects, proceeding from the lowest- 
mass main-sequence stars to lower effective temperatures is mainly characterised 
by the onset of condensate formation in the latest M dwarfs, leading to strongly 
dust-dominated atmospheres of early L-types with very red JHK-colours. From 
mid-L to early T- types this trend is reversed with a decreasing signature of dust, 
strengthening of water vapour absorption and the appearance of methane. The 
strongest CH4 bands at 3.3 and 2.2 urn have been observed in dwarfs as early as 
L5 • The L/T transition is thus now defined by the first appearance of the 
weaker band at 1.6 |im |7I8| . While the near-IR spectra of T dwarfs are turning 
blue again with decreasing T c ff, caused by the deepening molecular bands and 
collision-induced absorption (CIA) of H2 in the X-band, their optical-infrared 
colours become extremely red due to the massively pressure-broadened resonance 
lines of Nal and KI. 

Chemical interaction, the complex spectra of molecules, the physics of dust 
formation and the treatment of line broadening in high density conditions make 
spectral modelling of brown dwarfs and giant planets a challenging task. With 
advancements in computational techniques and availability of better chemical 
and physical input data over the years, models of the atmospheres and interiors, 
|13I14I15I16| have been able to keep up with observational progress in their 
ability to describe the global properties of brown dwarfs. In this review, we 
present an introduction into dwarf atmosphere modelling with the PHOENIX code 
in Sect. followed by an overview of the current status of models and recent 
developments in Sect. [31 We discuss the special case of modelling the irradiated 
atmospheres of close-in extrasolar giant planets in Sect.0]and discuss the outlook 
to possible direct observations of EGPs. 

2 The PHOENIX Code 

PHOENIX (see and references therein) is a general stellar model atmosphere 
code for treating both static and moving atmospheres, designed to be both gen- 
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eral enough to allow essentially all astrophysical objects to be modeled with a 
single code, and to make as few approximations as possible. Successful applica- 
tions of the PHOENIX code include models of Novae, all types of Supernovae ^Hj 
and Hot Stars ^Hj- Compared to these cases cool stars and brown dwarfs are 
in many respects simpler instances of the classical stellar atmosphere problem, 
i. e. the photospheres can be considered as plane parallel and the assumptions 
of hydrostatic equilibrium, local thermodynamical equilibrium (LTE) and flux 
conservation at each level hold to a high degree. Complications arise from the 
need to treat convection into optically thin layers, highly non-grey opacities that 
cause extreme deviations of the radiation field from that of a blackbody, and the 
prevalence of molecules at low temperatures. 



2.1 Atmospheric Structure and Chemistry 

Molecules first appear at temperatures of about 5 000 K, and dominate the condi- 
tions in M dwarfs, locking up hydrogen and most of the metals. This dramatically 
increases the complexity of the equation of state (EOS) as chemical interactions 
of hundreds of species have to be considered simultaneously. The PHOENIX code 
currently solves for the chemical equilibrium (hereafter CE) partial pressures 
of 40 elements with usually 2 to 6 ionization stages each, over 600 molecular 
species |14I20| . and over 1000 liquids and crystals from a study by Sharp & 
Huebner |5T]. Thermodynamic equilibrium models, however, can only describe 
the conditions, under which condensates may be formed, but not the actual nu- 
cleation rate and growth of grains. Therefore, two limiting cases for the presence 
of dust in the photosphere have been described by Allard et al. [23! , referred to 
as AMES-Dusty and AMES-Cond. While in both cases dust is assumed to have 
formed wherever it can under CE, removing refractory elements from the gas 
phase, the former models consider all grains to remain at the place where they 
formed, while the latter assume that dust will completely settle to deeper layers 
and does not contribute to the opacity at all. 



2.2 Opacities 

With the electronic, vibrational and rotational transitions of molecules being 
orders of magnitude more complex than atomic line spectra, the handling of up 
to several 10 8 line opacities poses a computational problem of its own. The input 
data size is illustrated by a breakdown of lines in our opacity database: 

• number of atomic/ionic spectral lines: 42 x 10 6 — * 0.6 GB 

• number of diatomic molecular lines (other than TiO): w 35 x 10 6 — ► 0.5 GB 

• number of TiO lines (1999): « 170 x 10 6 -> 2.5 GB 

• number of hot water vapour lines - before 1994: 0.035 x 10 6 



1994 
1997 
2001 



6 x 10 6 -> 0.1GB 
308 x 10 6 -> 4.5 GB 
100 x 10 6 -> 1.5 GB 



total molecular lines (Oct 2001): « 700 x 10 6 -> 10GB 
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Line blanketing is treated with the direct opacity sampling method, by dy- 
namically creating sub-lists of all transitions contributing to the opacity for the 
conditions of a given model, and calculating depth dependent Gauss (weak lines) 
or Voigt (strong lines) profiles [22] • For important resonance lines, such as the 
strongly broadened alkali Df/D2 lines, special damping profiles [22] are used. 
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Fig. 1. Combined dust extinction profiles for power-law grain size distributions with 1, 
2, 10 and 100 times the standard ISM values (full lines from bottom to top respectively, 
where the two first curves are nearly undistinguishable). The scattering and absorption 
contributions of the 100 ISM profile are also shown (dotted lines). The conditions are 
typical of the upper photospheric layers in a 1 800 K AMES-Dusty model atmosphere 
(see Allard et al. jZTjj for further details). The structures seen in the profile at A > 8.5/im 
are due to dust absorption (Mg2Si04 at 10 and 16.5 /im and MgA^CU at 13 /j,m). 
Scattering contributions dominate below 0.5/im, but remain modest in the IR for all 
but the largest grain sizes. The absorption profile shows little sensitivity to grain sizes. 

For models where dust opacity is considered (AMES-Dusty) , extinction from 
31 grain species are currently included, with absorption and scattering cross- 
sections calculated using the Mie formalism applied to spherical grains [21] ■ Fig- 
ureHshows the combined extinction from all dust species included in our EOS, 
which provide an essentially continuous opacity, though several pronounced fea- 
tures exist in the mid-lR, due to e. g. the absorption peaks of silicates and corun- 
dum (AI2O3). An uncertain factor in these calculations is the size of the grains, 
about which only little is known for the conditions of dense stellar atmospheres. 
For lack of a better model, a power-law distribution is frequently assumed, adopt- 
ing an ISM model of n(a) oc a~ 3 5 for grain diameters 6.25 nm< a <250nm. 
However, the total extinction per dust mass remains insensitive to particle size 
as long as the grains remain within the Rayleigh limit a <C A/2-7T, and therefore 
the IR opacity would only be affected if the grains became much larger than in 
the ISM (cf. Sect. IQl. 
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Fig. 2. Model spectra for the conditions typical of a VLM star (T e fr = 3 500 K, top, with 
very little dust being formed), a young BD (T c b — 1800K, middle, with dust opacity 
included) and an older field BD (T e a = 1000K, bottom, without dust opacity). The 
most important absorption features are indicated. 



3 Current Status of Models 

The evolution of spectral properties from the hydrogen-burning limit to the 
coolest observed brown dwarfs is displayed in FigJ^J The transition from M- 
to L- types is characterised by disappearance of the optical VO and TiO bands, 
while FeH becomes stronger in the red optical and near IR spectrum, and a shift 
of the IR SED, exhibited in increasingly red (J—K) colours. These changes are 
explained by depletion of gaseous TiO due to the formation of titanates and to 
a lesser degree the condensation of VO, and the added opacity of silicate grains. 
Dust extinction blocks efficiently the emission in the optical and tends to smear 
out the remaining molecular absorption bands, the most prominent being due 
to H 2 and CO. 

Dusty models such as the middle plot in Fig. |2l reproduce well the observed 
spectra of early L dwarfs. As T e g sinks below 1 700 K, silicates form at deeper 
layers, where they affect the emerging spectrum less, the outermost layers can 
cool and molecular bands begin to grow in strength again. The cooling also 
triggers a rapid change from a CE composition where most carbon is bound 
in carbon monoxide to a carbon chemistry dominated by methane, allowing 
the much stronger bands of this molecule to shape the IR spectrum. We can 
see the end point of this transformation at T e g = 1 000 K, with the bottom plot 
showing an AMES-Cond model where the atmosphere is dust-free, and water 
and methane bands are strong. Consequently, cool evolved brown dwarfs emit 
more flux at shorter wavelengths (in the near-IR Z, J, and H broad bands), 
between the strong water vapour and methane bands, than blackbodies of the 
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same bolometric luminosity would predict, and show very blue ( J — K) colours 
compared to L dwarfs. This important fact was already predicted by the model 
calculations of Allard et al. in 1995 ^3 an d confirmed by the discovery of the 
first evolved brown dwarf, Gliese 229 B |2l3j . 

3.1 Formation and Gravitational Settling of Dust Clouds 

Therefore, based upon observational evidence, the latest M and most L dwarfs 
carry atmospheric dust, whereas late T dwarfs show little or no evidence of con- 
densates contributing to the opacity. Our models presented in Sect . 12 . ll reproduce 
these extreme states well, but cannot adequately describe the transition between 
the two. In fact, in a static equilibrium model there is no solution that allows 
dust to remain suspended in the photosphere at all, as the grains will eventu- 
ally sink to lower layers, and take all condensible elements with them until their 
abundance has decreased below the saturation pressure. Radiation pressure is 
negligible at such low luminosities and certainly cannot provide the force re- 
quired to keep grains lifted up. Because of the relatively inefficient transport of 
energy by radiation and the large molecular opacities, these cool atmospheres 
are all convective. Convection can thus provide the essential replenishment of 
condensible material in the cloud layers. We therefore understand the presence 
of dust as a delicate equilibrium between 1) the sedimentation of the largest 
grains by gravity to deeper, hotter layers where they are melted and sublimated, 
2) the convective updraft of condensible material, and 3) the efficiency of the 
condensation process and growth rate of grains. It becomes clear from this that 
the cloud layer can only exist in close proximity to a convection zone, where 
mixing is still vigourous enough to balance the sedimentation and depletion pro- 
cess. As the radiative-convective boundary retreats with decreasing luminosity, 
the cloud top also sinks deeper into the atmosphere. 

Realising this connection between convective mixing and the presence of dust, 
various models have been developed to describe the cloud structure (cf. the 
review by Marley et al. (23)- The model of Tsuji |2H] assumes a fixed grain 
size and defines a lower limiting temperature T cr i t , below which grains would 
coagulate too fast to stay suspended. The vertical extent of the cloud deck is 
thus defined by the condensation temperature and T C rit, the latter being an 
empirical parameter that needs to be adjusted by comparison with observations. 

A more detailed approach has been developed by Ackerman & Marley \'27\ . 
which uses a simplified description of the equilibrium between grain settling 
and mixing to derive modal grain size and dust density. Above the convection 
zone, vertical transport is assumed to occur only due to turbulent mixing with 
a fixed eddy diffusion coefficient. Their model also includes a free parameter 
/ se d to describe the sedimentation efficiency, which again is adjusted to fit the 
cloud thickness to observations. They find that observed L dwarf colours can be 
reproduced well using a relatively small range in / se d |28| . 

In the model used in the current version of PHOENIX, Allard et al. |29IH()j 
describe a more detailed calculation of the equilibrium between mixing and var- 
ious condensate growth processes. Vertical transport within the convection zone 



Brown Dwarfs and Hot Jupiters 



7 



is calculated according to mixing length theory, and overshoot into the radiative 
layer is modelled by assuming a velocity law based on the 3D-hydrodynamical 
simulations of Ludwig et al. |31| . The time scales for condensation, coagulation 
and coalescence, and the sedimentation speed of grains, are currently calculated 
according to the description of Rossow |32| for cloud microphysics in planetary 
atmospheres. The removal of elements from the gas phase by condensation and 
sedimentation is self-consistently taken into account in the solution of the CE 
system. The main differences from |27] are thus the formulation of a depth- 
dependent mixing velocity driven by convective overshoot in the radiative zone, 
and the use of the microphysical time scale approach of Rossow, which adds 
a level of physical detail, though it risks extrapolating results which have only 
been tested empirically in the very different conditions of solar system planets. 
Another model based on the work of Rossow has been published by Cooper et 
al. |33| . though they give no details on their the mixing velocity law, and have 
not yet included the cloud opacity into the atmosphere structure. 

The most advanced and detailed treatment of cloud microphysics to date has 
been presented by Woitke & Helling 34, cf. also C. Helling, this volume]. Their 
model provides a kinetic description of turbulent mixing and dust formation and 
precipitation processes, also based on an overshooting velocity field calibrated on 
hydrodynamical models [HI]. At this stage only one species, TiO, is considered, 
and the computational complexity of these non-local cloud formation models has 
so far prohibited incorporation into a self-consistent model atmosphere code. 

It is worth noting that all these models except predict the formation of 
dust grains larger than lp.m, and ranging up to several 100 (Xm [33], in some 
layers of the atmosphere. IR extinction coefficients in these models are thus no 
longer independent of the grain size, and measurements of the near- and mid-IR 
SED may provide an observational test of the size distribution. 

In parallel to the settling into deeper layers it is also possible that the cloud 
deck is breaking up into smaller features, with the surface coverage decreasing 
with T e g. The balance of dust rainout and replenishment by convective upwelling 
is certainly conductive to the formation of unstable cloud structures, and evi- 
dence of this might be found in the very rapid change of the SED and the 
reappearance of FeH absorption bands in early T dwarfs [55| . But while obser- 
vations of both periodic and transient variability in L and some early T dwarfs 
clearly indicate the existence, and rapid evolution, of surface patterns 36 31], no 
trend for an increase in amplitude has been found that would support a signifi- 
cant change of the cloud fraction towards the L/T boundary [2E]. The models of 
Ackerman & Marley require a steady decrease in cloud coverage to explain the 
observed colour change at the beginning of the T dwarf sequence |25| , while the 
settling code of Allard et al. [2H| is able to reproduce these broad-band colours 
with a homogeneous cloud deck, but has not been tested against all changes of 
spectral features. More detailed comparison with spectroscopic observations will 
therefore be required to decide if uniform cloud models can fully account for the 
transition from dusty to dust-free atmospheres. 
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Fig. 3. Theoretical models of the strongest methane absorption band for a T dwarf 
of T e g = 1000 K. The plots, as labeled from the top, show spectra based on the 
GEISA/HITRAN line list, and on 3 lists calculated with the STDS code E2]. 

3.2 Molecular Abundances and Opacities 

A serious obstacle to improving models for ultracool atmospheres is the lack 
of physical input data for many atomic and molecular opacities. Most exist- 
ing databases of molecular lines such as HITRAN [35] have been compiled for 
planetary atmosphere studies, and are oriented towards low temperatures. Con- 
sequently, they are highly incomplete at energy levels strongly populated in 
brown dwarf and cool star atmospheres. Stellar atmosphere models therefore 
rely strongly on the efforts of laboratory spectroscopists and molecular physi- 
cists to produce comprehensive experimental and theoretical high-temperature 
line lists [241401 cf. also S. Johannson, this volume]. 

For dense molecular bands completeness of all relevant transitions (i. e. cov- 
ering all levels populated at temperatures of 500 - 3 000 K), is generally more 
important to correctly include blanketing effects, than accurate individual line 
strengths and positions. This is illustrated in Fig. [3| which compares lines from 
the low-temperature databases with simulations using the STDS code {41142) . 
calculated for transitions from the vibrational ground state, as well as the first 
(Dyad) and second (Pentad) excited levels (hot bands). For population num- 
bers typical of brown dwarf conditions, lines from each band contribute about 
the same integrated absorption strength, but distributed over approximately 
2 x 10 4 , 2 x 10 5 , 2 x 10 6 and 1.3 x 10 7 lines, respectively. As a result, the spectra 
calculated from the GEISA/HITRAN lines and from only the ground-state tran- 
sitions allow several times more flux to escape between the lines. The spectra 
including the hot bands on the other hand converge as soon as ~ 10 6 lines are 
included. 

Completeness and accuracy for the IR water bands in the coolest brown 
dwarfs remains problematic jllj . The available data for ammonia bands, which 
dominate the mid-IR spectrum in late T dwarfs, are even more limited. In addi- 
tion, mounting evidence suggests that CH4/CO as well as NH3/N2 do not reach 
their CE mixing ratios in the cool upper layers of T dwarf atmospheres, as some 
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Fig. 4. Models for a non-irradiated and irradiated (T e g = 500K, log g = 3.5) planet 
when located 1.0, 0.5 and 0.3 AU from a G2 primary star from I44| . Thermal structures 
for dust-free and dusty models are shown on the left, with the bottom line indicating the 
case without irradiation. The filled symbols refer to different optical depths (r) defined 
at A = 1.2 (Xm. The corresponding spectra are plotted on the right. For comparison, 
the spectrum of the G2V star is also shown. All fluxes have been scaled appropriately 
for the size of the planet and primary at a distance of 15 parsecs. 

key reaction times could be orders of magnitude larger than the timescale for 
mixing with warmer layers |43|. Probing molecular abundances at low spectral 
resolution is complicated by the strong blending of many bands (e. g. CO and 
H2O at 4.55u,m), which might be disentangled by high-resolution observations. 

4 Hot Jupiters 

The discovery rate of EGPs has only been slightly behind that of brown dwarfs. 
Most known EGPs are associated with a parent star similar to our Sun and 
have very small orbital separations compared to that of Jupiter. Their close 
proximity to a solar-type star makes the modeling of such objects slightly more 
complicated than brown dwarfs because both the intrinsic interior source of heat 
and the intense extrinsic solar radiation must be accounted for. Barman et al. 01] 
have modeled EGPs at various orbital distances and cooling stages. 

Figurc0]demonstrates the effects of irradiation on the atmosphere of an EGP 
which has cooled internally to 500 K. This corresponds to an age of about 10 
million years (Hj- This particular case illustrates the dramatic heating that oc- 
curs in the atmosphere of EGPs that are found as close as 0.05 AU from their 
parent star (hereafter "Hot Jupiters"). Stellar radiation is capable of warming 
the planet's surface to an equivalent temperature above 1600 K, i.e. conditions 
similar to dusty L dwarfs. Therefore, depending upon the atmospheric com- 
position history, we also may expect to find magnesium silicate clouds at the 
surface of "Hot Jupiters" . These models also demonstrate the effects of irradia- 
tion for the extreme cases of condensation in an EGP atmosphere. We compare 
the temperature structures for a cloud-free and cloudy planet located at various 
orbital separations. In the cloud- free model, the stellar radiation penetrates deep 
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Fig. 5. The luminosity ratio for a lMj EGP located 0.05 AU (top curve) and 5.0 AU 
(bottom curve) from a G2V main sequence star. In both cases, the EGP radius is 
1.42 _Rj and the intrinsic effective temperature was 500 K, typical of a young planet (a 
few Myrs) |5J. The location of various photometric bands are also shown. The contrast 
is most favourable beyond 4 urn. 

into the atmosphere and heats the layers near the radiative-convective bound- 
ary. At small orbital separations, the heating is sufficient to suppress convection 
in most of the atmosphere. In the cloudy case, the dust effectively shields the 
inner layers of the atmosphere and heating only occurs at the top of the atmo- 
sphere. The cloud-free EGP atmospheres have spectra similar to the dust-free, 
free-floating brown dwarf shown as the lowermost curve in Fig. However, as 
the planet is brought closer to the parent star, the molecular absorption features 
weaken, due to a flattening of the temperature profile. The dusty atmospheres 
have nearly featureless spectra due to the "grey" characteristics of dust opacity. 
When 0.05 AU is reached (not shown), the upper atmosphere is thermalised and 
the ncar-IR spectrum featureless in both cases. However, these models represent 
the substellar point (i. e. the point closest to the star) and therefore only predict 
a small portion of the complete observable spectrum. 

It is very likely that a large day/night temperature difference will exist lead- 
ing to a horizontal temperature gradient. Thus, the observable hemisphere of an 
EGP may have a complex spectrum composed of regions with varying degrees of 
stellar illumination and a fraction of the night side, depending on the phase and 
inclination. The exact thermal profile will also depend on the global atmospheric 
circulation that would redistribute the irradiated energy 05]. On the other hand, 
for orbital separations larger than about 1 AU the stellar flux no longer adds 
significantly to the thermal near-IR EGP spectrum, implying that such EGPs 
can be well approximated by a non-irradiated atmosphere. Such an approxima- 
tion is only valid for the thermal radiation emitted beyond 1 (J,m. At shorter 
wavelengths, the planet's atmosphere reflects the stellar flux, making it appear 
much brighter than isolated brown dwarfs. The reflected optical starlight per- 
sists out to large orbital distances, as cooler atmospheres form highly reflecting 
cloud layers (ices with albedos of 0.5). 



Brown Dwarfs and Hot Jupiters 



11 



Since planets have been most often found close to a star, it is important to 
identify the most favourable wavelengths for the direct detection of light from 
those planets. Figure presents the monochromatic luminosity ratio of an EGP 
at 0.05 and 5AU from a G2V star. Independent of the spectral type of the 
primary (from G to M), the EGP at 0.05 AU has a contrast reaching 10" 2 in the 
V band. However, when the planet is at 5 AU, the contrast reaches only 10 -3 5 at 
5|im, decreasing to ^10~ 5 in the H bandpass. Note that these numbers are, for 
reasons above, optimistic estimates and assume a "full moon" planetary phase. 
The contrast turns up at optical wavelengths due to starlight reflected by the 
planet. The night side of a "Hot Jupiter" may look more like the EGP with an 
orbital separation of 5 AU, except of course that there would be no reflection 
effects, with such a large day/night difference likely producing strong temporal 
variations of the contrast as the planet orbits the star. 

5 Conclusions 

Ten years of observational and theoretical work have brought us thorough under- 
standing of the basic physical properties of brown dwarfs. While many details, 
especially regarding non-equilibrium processes, and temporal and spatial vari- 
ability of the atmosphere, still need to be addressed, current models provide a 
solid foundation for understanding extrasolar giant planets, and predicting their 
observable signatures. BD and EGP (unirradiated as well as "Hot Jupiter" -like) 
model atmospheres, thermal profiles, synthetic spectra, photometry, and evolu- 
tion calculations are available for all evolutionary stages from a few Myrs from 
our web site at: http://perso.ens-lyon.fr/france.allard/ 
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